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Abstract 

We describe the resolvent approach for the rigorous study of the meso- 
scopic regime of Hermitian matrix spectra. We present results reflecting 
universal behaviour of the smoothed density of eigenvalue distribution of 
large random matrices. 

Random matrices of large dimensions introduced and studied by E.Wigner 
[ p0[ have applications in various fields of theoretical physics (see e.g. mono- 
graphs and reviews q, |], and references there). In these studies, the spectral 
properties of random matrix ensembles play an important role. 

Here the universality conjecture for large random matrices, formulated by 
F.Dyson ffj, fj}, is known as the most interesting and challenging problem. It 
concerns the asymptotically local spectral statistics, i.e. the functions that 
depend on certain number q of eigenvalues of random N x N matrix An and 
this number remains fixed when N — > oo. 

Loosely speaking, the universality conjecture states that the local statistics 
regarded in the limit N — > oo do not depend on the details of the probability 
distribution P(An) of the ensemble but are determined by the symmetries of the 
ensemble. For example, the expressions derived for local statistics of Hermitian 
ensembles are different from those of real symmetric matrices. 

Given Hermitian (or real symmetric) matrix An, the distribution of its eigen- 
values < . . . < is determined by the normalized eigenvalue counting 
function 



a N (X) = a(X;A N ) := #{A JVJ < \}N 
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or, equivalently, by the associated measure 

a N (A) = f PN {\) dA, A = (a, b) C R, 

J a 

with the formal density 

1 N 

MA) = ^$>(A-Af>). (1) 
i=i 

The function £i/v(A) is called the empirical eigenvalue distribution function. Re- 
garding (Tjv(Ajv), it turns to be the local spectral statistics when considered 
with the intervals of the length |Ajv| = 0(l/N) as N — > oo. 

In general, the local spectral regime is rather hard to analyse rigorously. The 
universality conjecture is supported mainly for those ensembles of random ma- 
trices that have explicit form of the joint probability distribution ttn (Ai , . . . , Ajv) 
of eigenvalues. Starting from ttn, the same expression for m-point correlation 
function is derived by Dyson for the circular ensemble of unitary random ma- 
trices (CUE) [§, by Mehta for CUE 0, by Pastur and Shcherbina for matrix 
models ensemble plj| . This expression is given by the determinant of m x m 
matrix with the entries {sin7r(ii — tj)/n(ti — tj)},i.j — l,...,m. The same ex- 
pression is derived in for a random matrix ensemble with the entries that are 
independent random variables, whose probability distribution is a convolution 
of the Gaussian distribution and the arbitrary one. 

Our principal goal is to examine the presence of universality of the spectral 
characteristics for those ensembles of random matrices, for which the explicit 
form of the joint eigenvalue distribution ttn is unknown. For example, random 
matrix with independent ±1 entries falls into this class. Our claim is that the 
eigenvalue density (1) smoothed over the intervals Ajv C R possesses the uni- 
versal properties as N — > oo provided the length ijy = | Ajv | satisfies conditions 
1 < l N < N. 

We determine the smoothing (or regularization) of (2) by the formula 

« A) -£ l + iv"(A-A-)^ "« A,)dV < 2 > 

and note that in this case 

R$(\) =Im Tr G JV (A + iiV- Q )^- 1 , 

where Gn{z) = (An ~ Z )~ X ■ 

According to the above definition, £j^(A) with a — 1 represents the asymp- 
totically local spectral statistics. The opposite asymptotic regime when a = 
is known as the global one. In this case the limit 

g{z) = lim Tr G N (z)N-\ |Im z\ > 
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if it exists, determines the limiting eigenvalue distribution cr(A) of the ensemble 
{A^v}; that is 



a(A) = lim a N {\), g(z) = / (A — z) dcr(A). 

Regarding the global regime, the resolvent approach developed in papers 
[ ft3[ Q is proved to be rather effective in studies of eigenvalue distribution of 
large random matrices (see, for example [l2|, |l^|). In this regime the limit of 
Qn{ z ) — Tr G^(z)N~ 1 depends on the probability distribution of the ensemble, 
i.e. is non-universal, as well as the fluctuations of gN( z ) f|, [T^ |. 

We are interested in the behaviour of (2) in the case of < a < 1. This 
regime is intermediate between the local and the global ones. It can be called 
the mesoscopic regime in random matrix spectra. For this regime, the modified 
version of the resolvent approach was proposed in to study spectral proper- 
ties of random matrices with independent arbitrary distributed entries (see also 
I !)• 



As the further development of the resolvent approach of 1 1 1 , we present the 
results concerning random matrices with statistically dependent entries. We 
consider the ensemble of random matrices 

^ m 

H m ,N(x,y) = —^2^{x)^(y), x,y = l,...,N, (3) 

where the random variables {^(x)}, x, fj, € N have joint Gaussian distribution 
with zero mathematical expectation and covariance 

(£)£„(?/)} = u 2 S xy 6^. 

Here S xy denotes the Kronecker delta-symbol. This ensemble first considered 
in is now of extensive use in the statistical mechanics of disordered spin 
systems jl8| and in the modelling of memory in the theory of neural networks 



Theorem 1. Let G m ^(z) — (H m jy — z) 1 . Then, for N,m — ^ oo, m/N - 
c > 0, the random variable 

^Sv(A) := Im Tr G rn , N {\ + iA^A^ 1 
converges with probability 1 as to the nonrandom limit 

1 



^c(A) - ^V4™ 4 -[A-(1 + CU 2 )] 2 (4) 
provided < a < 1 and A e A Cj „ = (u 2 (l - v / c) 2 ,u 2 (l + Vc) 2 )- 

Theorem 2. Consider k random variables 

7&(0 == RS(Ai) - E< ) Ar (A 4 )l ,i = l,...,k, 



3 



where Xi = A + TiN~ a with given ti. Then under the conditions of Theorem 1 
the joint distribution of the vector (7^(1), . . . ,7jv(fe)) converges to the Gaussian 
k- dimensional distribution with zero average and covariance 

Remark. It is easy to see that if \t± — t?\ — > 00, then 

C(tx,t 2 ) = -(r 1 -r 2 )- 2 (l + o(l)). (6) 

This coincides with the average value of the Dyson's 2-point correlation function 
for real symmetric matrices considered at large distances \t\ — ia| S> 1 0. 

To discuss these results, let us first note that theorem 1 proves existence of 
the smoothed density of eigenvalues that coincides with that derived in Jl6| in 
the global regime; g c (X) = a' c (X),X > 0, where 

cr c (A) = Jim a(X\H mN ). 
This density obviously differs from the semicircle (or Wigner) distribution 

<7 W (X) 

a w (X) = lim a(X;W N ) 



where Wn(x,v) = w(x,y)/yN are random symmetric matrices with indepen- 
dent identically distributed entries with zero mathematical expectation and vari- 
ance v 2 . This ensemble is known as the Wigner ensemble of random matrices. 
It is known since the pioneering work of Wigner |20| that 

m / m 1 / V4v 2 -A 2 , if|A|<2», 
Qw (X)=a w {X) = ^{ 0j ii\x\>2v: 

It should be noted that in papers j|, |ll| we have proved analogues of The- 
orems 1 and 2 for the Wigner ensemble of random matrices. We have shown 
that Theorem 2 is true provided Eiw(x,y) 8 is bounded and a £ (0,1/8). The 
correlation function C(ji,Tj) is given again by (5). Comparing these results, 
we conclude the fluctuations of the smoothed eigenvalue density do not feel the 
dependence between matrix elements. 

Thus, our results can be regarded as the statements corroborating the uni- 
versality conjecture for mesoscopic regime. Namely, they show that in the meso- 
scopic regime the smoothed density of eigenvalues is the selfaveraging vari- 
able. It converges as N — > 00 to the eigenvalue distribution of the ensemble and 
depends on the probability distribution of the random matrix ensemble. At the 
same time the fluctuations of in the limit N — > 00 coincide for two such 
different classes of random matrices as (3) and the Wigner one. 

This dual-type behaviour of the eigenvalue density of random matrices in 
the mesoscopic regime is well-known in theoretical physics (see, for example the 
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review ||, Chapter 8). For example, the universal properties of the mesoscopic 
eigenvalue density is studied in || for the matrix models ensemble. It was 
shown that the correlation function of the eigenvalue density (in theoretical 
physics terms, the "wide" correlator) depends on the edges of the spectrum. 
For the case of the symmetric support of the limiting eigenvalue distribution, 
the expression for the "wide" correlator coincides with the asymptotic expression 
(6). It should be noted that our result (6) does not depend on the support of 
QcW- 

Let us describe the method developed for the proof of Theorems 1 and 
2 (the full version will be published elsewhere). It represents a modification 
of the resolvent approach proposed in papers Jl3], [l4|. This approach was 
developed to study the eigenvalue distribution of random matrices and ran- 
dom operators in the global regime |Im z\ > as N — ► oo. It is based on 
derivation and asymptotic analysis of the system of relations for the moments 
= Ei[gpf(z)] k , k > 1. These relations are of the following form 

L k =aL k - 1 +bL k+l + <$>[ N \ (7) 

where the terms can be estimated by A _1 |Im z\~ k . This a priori estimate 
implies that in the study of the asymptotic behaviour of g^iz), one can restrict 
oneself with only two first relations. Namely, all information about the limiting 
behaviour of k > 1 can be derived from relations for L^' and 

To consider Qn(z) m the mesoscopic regime, we start with the same system 
of relations for ijj. • The main observation made in |]llf is that in this case 
we need all the infinite system of relations. More precisely, the more close a 
is to 1, the greater number K{a) is such that we need to consider relations for 

L{ N) ,k>K{a). 

The matter is that the terms can be estimated in terms of Lrj^ multi- 
plied by N" 13 , where (3 — min{a, 1 — a}. The structure of relations (7) is such 
that Lj, j < k enter into relation for L^' with the factor N@( k ~i\ There- 
fore, admitting a priori estimate |£x | < N a , we deduce that it enters into 
the relation (7) with the factor N~ k/3 . Regarding, k > K(a), one obtains the 
relations with the terms that converge to finite limits as N — > oo. 



Acknowledgments. The authors are grateful for A. Its, P.Bleher, and 
H.Widom and other organizers of the semester "Random Matrix Models and 
Their Applications" at MSRI (Berkeley) for the kind invitations to participate 
the workshops and for the financial support. 

References 

[1] J Ambj0rn, J Jurkiewicz, Yu M Makeenko, Phys. Lett. B 251 (1990) 517 



5 



[2] A Boutet de Monvel and A Khorunzhy, Markov Proc. Rel. Fields 4 (1998) 
175-197 

[3] A Boutet de Monvel and A Khorunzhy, Rand. Oper. Stoch. Equations 7 
(1999) 1-22 

[4] A Boutet de Monvel and A Khorunzhy, Rand. Oper. Stoch. Equations 7 
(1999), to appear 

[5] E Brczin and A Zee, Nucl. Phys. B 402 (1993) 613-627 

[6] F Dyson, J. Math. Phys. 3 (1962) 140-175 

[7] F Dyson; J. Math. Phys. 13 (1972) 90 

[8] T Guhr, A Muller-Groeling, and H A Wcidenmiiller Phys. Rep. 299 (1998) 

[9] J A Hertz, A Krogh, and R G Palmer. Introduction to the theory of neural 
computations, Addison- Wesley (1991) 

[10] K Johansson, preprint, (1998) 

[11] A Khorunzhy Rand. Oper. Stoch. Equations 5 (1997) 147-162 

[12] A Khorunzhy, B Khoruzhenko, and L Pastur J. Math. Phys. 37 (1996) 
5033-5060 

[13] A Khorunzhy and L Pastur, Commun. Math. Phys. 153 (1993) 605-646 

[14] A Khorunzhy and L Pastur, in: Spectral Operator Theory and Related 
Topics, Adv. Soviet Math. 19 (1994) 97-127 

[15] A Khorunzhy and G J Rodgers, Rep. Math. Phys., 43 (1998) 297-319 

[16] V Marchcnko and L Pastur, Math. USSR Sbornik 1 (1967) 457-483 

[17] M L Mehta. Random Matrices. Acad. Press, New York (1991) 

[18] M Mezard, G Parisi, and M Virasoro. Spin glass theory and beyond, World 
Scientific, Singapur (1987) 

[19] L Pastur and M Shcherbina, J. Stat. Phys. 86 (1997) 109-147 

[20] EWigner, Ann. Math. 62 (1955) 548-564 



6 



